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Abstract 

Hadronic matrix elements are expressed in the framework of BS formalism, and in 
the limit of heavy quark mass mg —> oo, we find that all form factors in processes 
of R —> D weak transition can be incorporated into a single universal Isgur-Wise 
{I — W) function Ru • v')^ which was first discovered by Isgur and Wise in heavy 
quark effective theory (HQET). By introducing a phenomenological interaction 
kernel, we solve the BS equation in leading order of l/mg, and calculate the 
numerical result for Ru • v'). We find that the slope of Ru • v') at the zero recoil 
point V ■ v' = 1 is larger for the process in which the involved light degree of 
freedom is s quark than the case of u,d quarks. Our results agree with others. 


In recent years a remarkable progress has been achieved in the physics of systems containing 
one heavy quark. HQET points out that in the limit of heavy quark mass mg oo a spin-heavy 
flavor symmetry [Q emerges and the weak decays of heavy hadrons can be described by a single 
form factor (the “/ — W function”) ' "^0 [H (where v and v' are the four-velocities of heavy 
meson before (after) the transition). Since I—W function plays a center role in the weak decays of 
heavy hadrons, many authors try to study its analytical properties approximately by htting with 
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the data. We list some of results as following: ref. 0 gives = 1.18±0.17; 

^ ^ — 

ref. @] gives = H-a(l —(u-u')^), a = 0.53; ref. p gives 5 = 0.91±0.03 

It can be easily derived from current conservation that ^{v ■ v') is normalized to unity at the 
zero recoil point. For B ^ D semileptonic decay processes the kinematic region is small due 
2 Mt,M Therefore if we know the slope at n ■ n' = 1, • v') can be 


, / /\ _ M?+M^ 

to (n ■ V )max 


approximately written as ^{v-v') = l+p{v-v' — 1) for small non-zero recoil. We can not predict the 
^{v-v') in theory only by symmetry consideration, because long-distance effection associated with 
strong interactions of light degrees of freedom so far can not be calculated from hrst principles. 
Therefore many models have been used to attempted to estimate the p. Lattice |0 and QCD sum 
rule have been applied to the calculation with rather different conclusions. Nonrelativistic 
quark dynamical models 0 appear to success in deriving the spectra of Qq systems but fail 
in discussing the form factor, they gave a p which is largely underestimated, because the light 
degrees of freedom can not be treated as light constituent quarks with nonrelativistic kinematics 
in Qq systems. Recently Close 0 has calculated p in a relativistic quark model, which agree 
with the empirical value for b ^ c transition. In his calculation, a gaussian ansatz wave function 
is used, they did not solve the bound state equation which associate with the interaction between 
quarks. 


In this paper we will use a relativistic quark dynamic model in which mesons are treated as 
a bound state consist of a pair of constituent quark and anti quark and described by relativistic 
BS wavefunction. In the limit of rriQ —*• 00, we hnd that S wave 0“ and 1“ states satisfy the 
same reduced BS equation. By introducing a phenomenological interquark potential we solve 
the zeroth order wave function at rest frame, which can be generalized to moving frame simply 
by the substitution q ^ q±. Then we apply this formalism to the calculation of form factors 
for weak transitions between heavy mesons and give the expressions of universal I — W function 
together with its slope at zero recoil point in terms of the solved zeroth order wave function. At 
last the numerical results are derived. 


We start from BS wavefunction of meson, which is dehned by 

Xp{.xi,X2) =< 0|T(Ti(a;i)T2(a;2))|i^ > • 


( 1 ) 


In the momentum space, it can be written as 
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( 2 ) 


Xp{q) = j d^qXp{x)e-^^-\ 

Here we use the standard center of mass and relative variables 

X = T]iXi + r]2X2, X = Xi - X2, (3) 

where rji = {i = 1,2). Then the bound state BS equation in momentum space reads 

{:^i-mi)xp{q ){^2 + m 2 ) = ^ J d*kG{P,q - k)xp{k), (4) 

where pi and p 2 represent the momenta of quark and antiquark repectively, G{P, q — k) is the 
interaction kernel which dominates the inter-quark dynamics. According to Eq.(^ we have 

Pi = 7]iP + q, p2 = P 2 P - q- (5) 


We define the three dimension BS wavefunction 


*F® = / MdXpii), 


( 6 ) 


which satishes the BS equation iTl[| : 

4>f(9) = ^ 


po-El- Ep 


/ d^kG{P,q-k)^p{k)YAi, 


( 7 ) 


where G{P, q — k) represents the “instantaneous” part of the potential G(P, q — k), and and 
A?, are positive energy projectors of the quark and antiquark respectively, 

{El + 7o7 ■ Pi + mixo) 


Al = 


2Ei 


( 8 ) 


^2 ^ (^2 - 7o7 ■ P 2 - m27o) 

Here Ei = \Jrn\+ pi, E 2 = \jrn 2 P Pi- deriving the equation (0) we have used the “instanta¬ 
neous approximation” and neglected the negative energy projectors in the quark propagators. 

In the moving frame of velocity v, the BS wavefunctions of S wave 0“ and 1“ mesons can be 
written as following: 
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= A+7o(1 + ^)757oA^(/?(P, g); 


( 10 ) 


= A+7o(1 + 


( 11 ) 


where = ^, and M are the 4-momentum and mass of meson. While ^ = 7 ^e^, where 
is the polarization vector of 1“ meson. (p(P, ^ and /(P, g) are scalar functions of variables 
P and q in general. Here for convenience we assume that ip and / only depend on variable 
|g±| = |g ~ gill = |g ~ (g • 'y)'y|. This approach is equivalence to decomposing relative momentum 
variable g into gy and g^ instead of go and g in the moving frame. It is obviously that in the rest 
frame v = (1,0), |g±| = |^. 

In the heavy quark limit mi —>■ oo, we get Pi —>■ P^, p 2 = m 2 V — g, and 

A^Pi)7o = (12) 

ZVq 


7oA^(P2) 


p27o -^■p2-m2 

2 E 2 


Thus the wave functions ([I0| ) and (|^ read 

= ^!^A757oA7(I?±I); 


(13) 


(14) 


Kf® = (15) 

^ 2no 

The interaction kernel is usually taken as the one consisting of a scalar and a vector kernel, which 
can be written as 


G(g-)$(g-) = Gsmi^ + Gv{q)7,mr 


(16) 


Substituting ([T^), (HH) and ([l6|) into (|^) and making v = (1,0), we obtain that ip{\^) and /(|^) 
satisfy the same reduced BS function. 


60o-(k1) = ^20o-(|g1) - / d^k{Gs{q- k) - Gv{q- k)) 


E 2 + m 2 
2 E 2 


(l>o-{\k\) 


f £‘k(Gs(q- k) + Gv(q- k))k:^Mklh 

tj2 + m2 J zrlj2 


(17) 
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where e = M — mi. We can easily obtain the scalar wave fnnction (p{\q±\) and f{\q±\) in moving 
frame only by making variable snbstitntion q q± in (^(|(?|). 


In this BS description, the hadronic matrix elements of weak transition \P,M >—> |P', M' > 
are given as El 


< P', M’\Q'TQ\P, M >= / d'p 2 TrxA']')rxp{a){h + m 2 ), 


(18) 


where Xp'W) = ~7°Xp'(5'07°- Eq.fpISD can be rednced into a simpler form which is expressed 
in terms of the three dimensional BS wave fnnctions, on condition that the negative energy 
projectors in qnark propagators are neglected and the kernel is independent of the relative time 
variable 

< P\M'\QTQ\P,M>= {27r)^ I (19) 


where the qnark may change its flavor while the antiqnark remains a spectator. Then the 
normalization of the wave fnnction reads 

{27rf j = 2E = 2\J+ P^ (20) 


It is well known that, in the heavy qnark effective theory, we have the following relations for 
B D weak transition in the limit mg —> cx): 

(i). Matrix element of vector cnrrent for 0“ —0“, 

< p'\v,\p >= vmpav • v'Kv, + vu-, ( 21 ) 


(ii). Matrix element of vector cnrrent for 0 —> 1 , 

< P'\V^\P >= -tVMAP^v ■ n')£^.a/3eV“n^; (22) 


(hi). Matrix element of axial vector cnrrent for 0 ^ 1 , 

< P'\A^\P >= V MM'^{v ■ n')[(l + V ■ v'Cfp — e ■ vv'^. (23) 

Here v and v' are fonr velocities of mesons before and after the transition. We hnd that all 
form factors are rednced to a single I-W fnnction ^{y-v'), which satishes normalization condition 
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^(1) = 1. I-W function is relative to the light degree of freedom in heavy meson and can not be 
determined only by symmetric analysis. In our approach we will easily derive the unitary I-W 
function, which can be expressed in terms of the leading order BS wavefunction. Additionally by 
introducing a phenomenological interaction kernel, we solve the BS wavefunction and calculate 
I-W function. 

Substituting the expressions of wavefunction (^) (|^ into the transition formula ( 0 ) , we 
hnd 

(i) . Matrix element of vector current for 0“ —> 0“, 

< -P'lnl^’ >= (2>r)“/<iy2rr|<*5; + (gi)7V<*t (,-)] = A(v^ + vW, (24) 

(ii) . Matrix element of vector current for 0“ —> 1“, 

< P'\V^\P >= {27rf I (^)]= (25) 


(hi). Matrix element of axial vector current for 0 ^ 1 , 

< P'\A^\P >= ( 27 r)^ [ d^p 2 Tr[^^p,^{^)'y°'y 5 'yf,^°p (^]= A[(l -F v ■ v')e^ - e ■ vv' ], ( 26 ) 


where 


A = 


2(27rf 

vOy/O 


m2 




E 2 1 + V ■ v' 


[nO -f n'O - pAv+v) 

L £/2 


]M\q±\MW±\) 


They are consistent with the result derived from HQET, comparing with (PT|), 
get 


^(n-n') = 


A 


2(27r 


VMAP VMWvov'o 

0 , ,0 P2-(v + P) 


'tJn'Hn J ^ 9 . 


-t- n'' 


E 9 


2 I + v -v' 


The normalization conditions of initial and hnal wavefunctions read 


f ^3 m2 + E2V^-p2-V ^ ^M^ + P^ 

{2tiY / Pp 2 --- (P {q±) = 

J V 0 E 2 


and 


(27) 


we 


(28) 


(29) 
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{27lfJ d^P2- 


^3 f ^ 3 ^ ^2 + E2v'^ - 4 :^^ 


M'2 + P'2 


v',E2 ^ 2 

For convenience, we assume initial state is rest v = (1,0),n' = (no,'F),then 


(30) 


I /O 

V ■ V = V = 


M'2 + p>2 


M '2 


I-W function can be simplified as 

2(27r)^ f ,3 ,E2 + m2 




( 31 ) 


Eo 


Now we remark about the relations of momentum variables appearing above. The usual 
quark model assume that the four momenta of the light degrees of freedom are conserved during 
the transition implies that p 2 = p '2 and hence that m 2 V — q = m 2 v' — g', so we get 

ql = q^-{q-vf = -(f = -if2, (32) 


q'l = - (g' • v'Y = {ni 2 v' - P 2 Y - {{m 2 v' - P 2 ) ■ v'f. (33) 


In order to get the relation of q\ and g(^ in terms of m 2 ,v,v', we must determine the zero com¬ 
ponent P 2 . If we go to the zero-binding limit, the light quark is on-shell P 2 = P 2 = (.^ 2 ,^ 2 ), E 2 = 
\Jm 2 + pi- Then 

= P 2 - (^2^0 -P2-V) . 


We employ the following interquark potentials including a long-ranged conhnement potential 
(Lorentz scalar) and a short-ranged one-gluon exchange potential (Lorentz vector) 


Vs{r) = \r 


1 - e" 


ar 


VAr) = 

3 r 


(34) 


where the introduction of the factor is to regulate the infrared divergence and also to incor¬ 
porate the color screening effects of the dynamical light quark pair on the QQ linear confinement 
potential. In momentum space the potentials become 


Gsip) 


--S^{p) + 

a 


X 1 

^2 (jp. _|_ q , 2^2 ’ 
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Gvip) = 


2 as{p) 


(35) 


3^2 _ j _ q ,2 ’ 

where as(p) is the quark-gluon running coupling constant and is assumed to become a constant 
of 0(1) as —>• 0 

127r 1 


(^sip) = 


27 ln{a + p‘^/AlcD)' 


(36) 


which saturation value is 


as{p = 0 ) = 


127r 1 


27 Ina 

The constant A, a, a and Aqqd are the parameters that characterize the potentials. Based on 
the kernel above we can solve Eq.(|T7|) numerical. We use 


A = 0.20Gen^, a = 1.0, Aqcd = 0.21Gen, a = O.OQGev, 


(37) 


and quark mass 


= 0.35Gen, m<j = 0.5Gen. 


With these parameters the correct 2S — IS* spacing and s — u,d, 1“ — 0“ splitting can be derived. 
Substituting (p{\q±\) and obtained from solving Eq. ([T7D into (^l|) , we get I-W function 

^{v ■ v') in the range of 1-1.6 for v ■ v' (Fig.l). We hnd they depend on the light degrees of 
freedom, the slope decrease as the mass of light quark increase, this is agree with the results of 
others 0. 


Finally we give the slope of I-W function at zero recoiling by differential of ^{v-v') for v-v' = 1 


P = 


di{v-v' )^ 
d{v ■ v' 
-2 


\w'=l 


d P 2 [ 


2(^2+ m 2 ),,,, :,^2-n.2 , P>^i\P2\) 


Eo 


-{i^\u;)r{p2-v'yEi + 


+P{\P2\)y'{^){P2-v')% 


(38) 


where 

/• 3 E 2 + m 2 
C = J d^P2 - — - (p (|P2|)- 

Here we have dehned a new function -0 satishes 

V’(^) = T'dgnI) = p{\P2\), 
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for a; = Using the wavefunction solving from ([T7|) with parameters (|^) , we get 


p„,rf = -1.14, p, = -1.22. 

In order to examine the sensitivity of calculated result to the potential parameters, we choice 
another group of parameters 


A = 0.2, a = 1.0, Aqcd = 0.15, a = 0.06, mu,d = 0.35, = 0.5 

we then get 

p„,rf = -1.18, p, = -1.27. 

We hnd that the value of p increase as Aqcd increase, although the change is very small. This 
is resonable because p mainly depend on the overlap integral between wavefunctions before and 
after transition <p(|q'±|)<p(|q'a|), when light quark mass increase or Aqcd decrease the wavefunction 
in momentum space become “thinner” and fall off faster, so supress the overlap part and enlarge 
the value of p. Our approach, in obtaining the scalar wavefunction at moving frame only by a 
convariant transformation ^ ^ in the wave function at rest frame, is resonable. 

In conclusion, our approach not only give the model dependent results in agreement with 
HQET in the heavy quark mass tuq oo, but also develop a model to calculate the I-W function 
which can not be predicted in HQET. We obtain the slope p = —1.1, p = —1.2 for light quark 
being u,d and s respectively, which is consistent with empirical value and have signihcantly 
changing comparing with the result of nonrelativistic model. The numerical values for p are 
insensitive to the choice of parameters in potentials. These show that in Qq systems, the light 
constituent quark is relativistic particle and relativistic effects must be taken into consideration. 
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Figl. Isgur-Wise function in weak transition process B ^ D with light quark u, d (solid) and 
s (dashed) 
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